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We investigate the potential to use a magneto-thermo-electric instability that may be induced in 
a mesoscopic magnetic multi-layer (F/f/F) to create and control magnetic superstructures. In the 
studied multilayer two strongly ferromagnetic layers (F) are coupled through a weakly ferromagnetic 
spacer (f ) by an "exchange spring" with a temperature dependent "spring constant" that can be varied 
by Joule heating caused by an electrical dc current. We show that in the current-in-plane (CIP) 
configuration a distribution of the magnetization, which is homogeneous in the direction of the 
current flow, is unstable in the presence of an external magnetic field if the length L of the sample 
in this direction exceeds some critical value L c ~ 10 fim. This spatial instability results in the 
spontaneous formation of a moving domain of magnetization directions, the length of which can be 
controlled by the bias voltage in the limit L > L c . Furthermore, we show that in such a situation 
the current-voltage characteristics has a plateau with hysteresis loops at its ends and demonstrate 
that if biased in the plateau region the studied device functions as an exponentially precise current 
stabilizer. 

PACS numbers: 



I. INTRODUCTION 



A useful tool for manipulating the local magnetic or- 
der in artificially structured materials is provided by the 
possibility to control the magnetization of a nanomag- 
net by injecting an electrical current. Several such sce- 
narios have been discussed, including the so-called spin 
torque transfer (STT) technique based on the suggestion 
by Slonczewski 1 and Berger 2 to use an injection current 
of spin-polarized electrons. The high current densities 
needed in this case can easily be achieved in electrical 
point contacts of submicron size, where densities of the 
order 10 8 -r- 10 10 A/cm 2 can be reached without signif- 
icant heating of the material,'*' 4 but for larger contacts 
thermal heating can not be avoided. Instead, Joule heat- 
ing caused by an (unpolarized) current can be used for 
thermal manipulation of magnetization. In Refs. 5 and 
6, this approach was proposed for varying the strength of 
the exchange coupling of two strongly ferromagnetic lay- 
ers (F) separated by a weakly ferromagnetic spacer (f). In 
Ref. 6, the ability of an external magnetic field to change 
the relative orientation of the magnetization in the outer 
layers of an F/f/F tri-layer magnetic stack was demon- 
strated. The result was that by varying the temperature, 
and hence varying the strength of the exchange-spring 
coupling through the spacer layer f, the relative orien- 
tation of the magnetization of the outer F-layers could 
be continuously and reversibly changed from being par- 
allel to being antiparallel. Consequently, Joule heating 
by forcing a dc current through the structure allows an 
electro-thermal manipulation of the relative magnetiza- 
tion directions. 



This kind of dc current-induced manipulation of the 
magnetization direction was further studied and observed 
in stacks predicted to have nonlinear N- and S-shapcd 
current-voltage characteristics (CVC) in Refs. 6 and 7, 
where temporal oscillations of the magnetization direc- 
tion, temperature and electric current in the magnetic 
stack were also investigated. 

In this work we explore the possibility to use Joule 
heating by a current-in-plane (CIP) dc electrical current 
to control the spatial distribution of the magnetization 
directions in an exchange-spring layered structure of the 
type 8 sketched in Fig. 1. We will show that, if the volt- 
age bias exceeds a critical value for which the sign of 
the differential resistance becomes negative, two coupled 
magnetic domain walls can spontaneously appear along 
the current flow (see Fig. 6) at a distance from each other 
that can be controlled by the bias voltage. 

Domain formation is of course not a new phenomenon. 
The Gunn effect, 9 well known from the physics of semi- 
conductors, is the name given to the spontaneous forma- 
tion of (moving) electric domains in a semiconductor bi- 
ased in a region of negative differential resistance (which 
requires an N-shaped CVC). Electric domains in normal 
metals can also appear and may be due to structural 1 " 
and magnetic 11 ' 12 transitions, a sharp temperature de- 
pendence of the resistance at low temperatures and mag- 
netic breakdown, 1,3-1 ' evaporation 1 '" - and melting 19 (for 
a review see, e.g., Ref. 20). In all these cases, however, 
the domain sizes are macroscopically large, typically sev- 
eral cm. In contrast, we will show here that the N- and 
S-shaped CVC:s of the magnetic exchange-spring struc- 
tures suggested in Refs. 5-7 may give rise to magneto- 
thermo-electric domains with a characteristic size of the 



2 



order of 10 fim. The spatial distribution of the magne- 
tization in these stacks can vary from one corresponding 
to a single magnetic domain wall to a spatially periodic 
magnetization distribution. 

The structure of the paper is as follows. In Sect. II we 
briefly discuss some important features of the tempera- 
ture dependence of the magnetization orientation in the 
exchange-coupled stack sketched in Fig. 1 and derive the 
N-shaped CVC that is the prerequisite for the magneto- 
thermo-electric instability discussed in Sect. III. There we 
show that under certain conditions an instability leads to 
spatially highly inhomogeneous distributions of the mag- 
netization direction in one of the layers (layer 2 in Fig. 1), 
of the temperature, and of the electric field inside the 
magnetic stack, corresponding to the formation of a sta- 
ble magneto-thermo-electric domain (MTED) structure 
in the stack. In the concluding Sect. IV we summarize 
the main results of the paper and estimate the values of 
the relevant parameters that lead to MTED formation. 
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FIG. 1: Orientation of the magnetic moments in the stack of 
three ferromagnetic layers discussed in the text. The magnetic 
moments in layers 0, 1, and 2 (whose relative magnitudes are 
indicated by the thickness of the short arrows) are coupled 
by the exchange interaction thus forming an exchange spring 
tri-layer. An external magnetic field H is directed antiparallel 
to the magnetization in layer and a current J flows in the 
plane of the layers (along the x-axis). 



II. N-SHAPED CURRENT- VOLTAGE 
CHARACTERISTICS OF A MAGNETIC STACK 
UNDER JOULE HEATING 

The magnetic stack under consideration has three fer- 
romagnetic layers as shown in Fig. 1. The outer two 
layers (0 and 2) are strongly ferromagnetic and coupled 
via the exchange interaction through a weakly ferromag- 
netic spacer layer (1). The Curie temperature of 
layer 1 is assumed to be lower than the Curie tempera- 
tures Tc°' 2 ^ of layers and 2. In addition we assume the 
magnetization direction of layer to be fixed. A static 
external magnetic field H, directed opposite to the mag- 
netization of layer 0, is required to be weak enough that 
at low temperatures T the magnetization of layer 2 is 
kept parallel to the magnetization of layer due to the 
exchange interaction between them via layer 1 . At H = 
and T > Tc this tri-layer is similar to the spin-flip "free 
layer" widely used in memory device applications. 21 The 
stack is incorporated into an external circuit in such a 
way that a current J flows through the cross-section of 
the layers and 
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Here R(Q) and i?o are the magnetoresistance and the 
angle-independent resistance of the stack, is the angle 
between the magnetization directions of layers and 2, 
and V is the voltage drop across the stack. 

In Ref. (> it was shown that a magnetic configuration 
with parallel orientations of the magnetization in layers 
0, 1 and 2 becomes unstable if the temperature exceeds 
some critical temperature T c (or) < T c (1) . The magnetiza- 
tion direction in layer 2 smoothly tilts with an increase 
of the stack temperature T in the temperature interval 
T (ot) < T < r a) The dependence of the equilibrium tilt 



angle between the magnetization directions of layers 
and 2 on T and the magnetic field H is determined by 
the equation 6 



= D(H,T) sinQ, T < T c (1) 
= ±tt, T > T c (1) , 
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Here L12 and M\^{T) are the widths and the magnetic 
moments of layers 1 and 2, respectively; a,\ ~ Ji/aAf 2 (0) 
is the exchange constant, I\ is the exchange energy in 
layer 1, fi B is the Bohr magneton, k B is Boltzmann's 
constant, and a is the lattice spacing. D(H,T) is a di- 
mensionless parameter that determines how effective the 
external magnetic field is to cause the misorientation ef- 
fect under consideration. More precisely, it is the ratio 
between the energy of magnetic layer 2 in the external 
magnetic field and the energy of the indirect exchange 
between layers and 2 (see Fig. 2). At low tempera- 
tures the indirect exchange energy prevails, the parame- 
ter D(H 7 T) < 1 and Eq. (2) has only one root, = 0, 
thus a parallel orientation of magnetic moments in lay- 
ers 0, 1 and 2 of the stack is thermodynamically stable. 
However, at temperature Tj or ' < Tc , for which 

D{T^ m \H) = 1, 

two new solutions, = ±|# min | 7^ 0, appear. The par- 
allel magnetization corresponding to = is now un- 
stable, and the direction of the magnetization in region 
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FIG. 2: Normalized temperature dependence of the mag- 
netic stack resistance, R — R[@(T)], as determined by the 
T-dependence according to Eq. (2) of the angle Q(T) be- 
tween the magnetization directions in layers and 2 of Fig. 1. 
The calculations were made for R~/R+ = 0.2, Do = 0.2; 
R± = R(tt) ± -R(O) and T c is the Curie temperature of the 
spacer layer 1 between layers and 2. 



2 tilts with an increase of temperature in the interval 
Tc (or) < T < T «_ Thc critical temperature T c (or) of this 
orientational phase transition is obtained from Eq. (3) 

as 22 



T (or) = T (l) j 



-.(1) 



D (ff). 



(5) 



The orientational transition discussed above can be 
detected by measuring the temperature dependence of 
the stack magnetoresistance, 21 R = R[Q(T)], plotted for 
a typical case in Fig. 6. This temperature dependence 
is caused by the temperature dependence of the mis- 
alignment angle O = 0(T), which is implicitly given by 
Eq. (2). 

If thc stack is Joule heated by a current J its temper- 
ature T(V) is determined by the heat-balance condition 



JV = Q(T), J = V/R eS (Q), 



where 



R(Q)Ro 



(6) 



(7) 



R(G) + Ro ' 

in conjunction with Eq. (2), which determines the tem- 
perature dependence of 0[T(V)]. Here V is the voltage 
drop across the stack, Q{T) is the heat flux flowing from 
the stack and i? e ff(0) is the total stack magnetoresis- 
tance. Here and below we neglect the explicit dependence 
of the magnetoresistance on T since we consider a thin 
stack in which elastic scattering of electrons is the main 
mechanism of the stack resistance. On the other hand, 
we consider the temperature changes caused by the Joule 
heating only in a narrow vicinity of Tc X \ which is suffi- 
ciently lower than both thc critical temperatures Tj ' 2 ^ 
and the Debye temperature. 

Equations (6) and (2) define the CVC of the stack 

w = R^wyyy < 8 > 



where Q(V) = Q[T(V)]. Below we will simplify the no- 
tation by dropping thc subscript "eff" and simply write 
R(G). 

From Eq. (6) it is clear that the differential conduc- 
tance of the stack, dJ/dV, is negative if the sample is 
Joule heated to a temperature at which 

d{Q(T)/R[@(T)]}/dT < 0. 

Therefore, the main properties of the system under con- 
sideration are determined by thc behavior of thc function 
X = Q(T)/R[Q(T)] plotted for typical parameters in the 
top panel of Fig. 4. Using Eq. (2) one may express the 
differential conductance in terms of the dependence of 
the magnetoresistance R on the magnetization angle 
as 6 



dJ ° _ fffe ^ [^~ 1 ( e )( 1 -- 6sine / e )] / 
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where [. . .]' means the derivative of the bracketed quan- 
tity with respect to 0, and 
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It follows from Eq. (9) the differential conductance 
dJ/dV is negative if 
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In this case the current-voltage characteristics (CVC) of 
the stack is N-shaped as shown in Fig. 3. 

Using Eqs. (1) and (9) and assuming the magnetore- 
sistance to be of the form 2 ' i?(0) = R+ (1 — rcos0), 
where 

r=^> 0; R±= m±m 

R+ ' 2 



(11) 



one finds that thc differential conductance dJ/dV is neg- 
ative if 



3r 
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R(0) - (1 - r) 2 R A 



(1 + 2r)R(0) + (1 - r) 2 (l - 4r)fl H 



(12) 



It follows that a CVC with a negative differential resis- 
tance is possible if R(0) > (1 — r) 2 R + . 

In the case that the stack has a negative differential re- 
sistance, nonlinear current and magnetization-direction 
oscillations may spontaneously arise if the stack is incor- 
porated in a voltage biased electrical circuit in series with 
an inductor. (> '' In this paper we show that another type 
of magneto-electrical instability can arise in such a stack 
if the electrical current flows in the plane of the layers 
(CIP-configuration): a homogeneous distribution of mag- 
netization direction, temperature and electric field along 
the spring-type magnetic stack becomes unstable and a 
magneto-thermo-electric domain spontaneous arise in the 
stack. Here and below we consider the case that the elec- 
trical current flowing through thc sample is lower than 
the torque critical current and hence the torque effect is 
absent. 
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FIG. 3: Current-voltage characteristics (CVC) of the mag- 
netic stack of Fig. 1 in which the magnetization directions 
in the layers are homogeneously distributed along the x- 
direction (along the stack). It was calculated for R(Q) — 
R+ - i?_cos9, R-/R+ = 0.2, D = 0.2; J c = V c /R{n); 

V c = \/i?(7r)Q(ri 1) )/^t- The branches - a and b - b' 
of the CVC correspond to parallel and antiparallel orienta- 
tions of the stack magnetization, respectively (the parts a — a' 
and — b are unstable); the branch a — b corresponds to 

o < e[T(v)] < 7T. 



III. MAGNETO-ELECTRO-THERMAL 
INSTABILITY IN A MAGNETIC STACK 

In this section we will work in the voltage bias regime, 
where the resistance of the external circuit into which the 
magnetic stack is incorporated can be neglected in com- 
parison with that of the stack. In this case, using the 
known relation between the electric field and the tem- 
perature (see, e.g., Ref. 28) and taking into account that 
the temperature, T(x,t), being a function of the coordi- 
nate x along the stack and time t, satisfies the continuity 
equation for the heat flow, one obtains a set of basic 
equations for the problem, 
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f(T,j) 
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da dT d 
df~dx ~ dx 
Q(T)/n st -j 2 (t)p{6(T) 
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-f(T,j) 



(13) 



where the T-dependence of 0(T) is given by Eq. (2). Fur- 
thermore, j(t) = J(t)L/S is the current density, which is 
independent of x due to the condition of local electrical 
neutrality, 5* is the cross-section area of the stack and L 
is its length, c v is the heat capacity per unit volume, a is 
the proportionality coefficient between the electric field 
and the temperature gradient," 8 k is the thermal conduc- 
tivity, Q(T) is the heat flux flowing from the stack, il st 
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FIG. 4: Plots of the function X (T) = Q(T)/R[G(T)] and 
the "potential" energy W(T, J) calculated for R(Q) = R+ — 
i?_cos6, R-/R+ = 0.2, D = 0.2, and J/J c = 1.03; J c = 
V c /R{tt). 



is its volume and p[9] = R e g[9]S/L is the stack magneto- 
resistivity [see Eq. (7)], the brackets (...) that appear in 
the last part of Eq. (13) indicate an average over x along 
the whole stack of length L. 

The boundary condition needed to solve Eq. (13) is 
the continuity of the heat flux at both ends of the stack 
(which is coupled to an external circuit with a fixed volt- 
age drop V over the stack) . We shall not write down any 
explicit expression for this condition, since it will become 
clear below that the magneto-thermal domain structure 
does not depend on the boundary conditions if L is suf- 
ficiently large. Instead, for the sake of simplicity, we use 
the periodic boundary condition T(x + L,t) = T(x,t). 

The set of equations (13) always has the steady-state 
homogeneous solution 



t = t (v), e = e„(v) = e[T (y)], 

j = p- 1 [Q(T )}V/L, f(T o ,j o ) = 0. 



(14) 



The last equation in (14) is identical to the energy bal- 
ance condition (6) which, together with Eq. (2), deter- 
mines the N-shaped CVC shown in Fig. 3. 

As shown in Appendix A, if the differential conduc- 
tance dJ/dV is negative the uniform magnetization along 
the stack is stable only if the stack length L is shorter 
than some critical length L c , where 



L c = 



p- 1 (e )(sine/e)' 



p- x (e) (i - L>sme/e) 



(15) 
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FIG. 5: Current-voltage characteristics (CVC) of the mag- 
netic stack of Fig. 1 for L ^> L c , i.e. when the stack con- 
tains a magneto-thermo-electric domain. The CVC was cal- 
culated for R(Q) =R+-R- cos 6, R-/R+ = 0.2, Do = 0.2 

V c /R(tt); V c = 



Jc 



a and b 



7?(7r)Q(T c (1) )/n st . The branches 
b' of the IVC correspond to parallel and an- 



tiparallel orientations of the stack magnetization, respectively, 
while the solid branch a — b corresponds to a stable magneto- 
thermo-electric domain (MTED) spontaneously formed inside 
the stack. The dashed lines indicate sections where the CVC 
is unstable. The red vertical arrows indicate the hysteresis 
loop in the CVC caused by the disappearance and emergence 
of the MTED as the bias voltage is changed. The thin hori- 
zontal line shows the stabilization current Jo. 



and the derivatives with respect to are evaluated for 
= Go- However, when L > L c the uniform distribu- 
tion of the stack magnetization becomes unstable against 
fluctuations comprising an arbitrary sum of harmonics 
A n ex.p(i2nnx / L) (n = ±1,±2...) with |n| < L/L c . A 
fluctuation with |n| > L/ L c or |n = 0| (uniform fluctua- 
tion), on the other hand, does not destroy the stability of 
the homogeneous solution (14) [see Eq. (A2)]. We note 
here that the characteristic value of the critical length 
can be rather short. Using Eqs. (15) and (9) and the 
Lorentz ratio n/a = tt k%T /3e, where a = 1/p and e is 
the electron charge, one finds that 
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for a realistic experimental situation: 



0.1 4- 0.3, 



(i) 



100 K, p(0) ~ 10/xQcm, j ~ 10 6 -j- 10 7 A/cm 2 



Therefore, in the range of parameters L > L c homo- 
geneous distributions of the magnetization direction, 0, 
temperature 



sin© 




[see Eq. (2) and Eq. (3)], and the electric field, 



£ = p[Q(T)]j 



(17) 



(18) 



along the system are unstable, and a magneto-thermo- 
electric domain, moving with a constant velocity s may 
spontaneous arise inside the magnetic stack: 

Q(x,t) = Q d {x - st) = Q[v(x - st,j d )], 

T(x,t) = v(x - st,j d ), j d (V) = (19) 

where the definition of the brackets (...) is the same as 
in Eq. (13), and v{x) satisfies the equation of motion of 
a fictitious particle of "mass" k sa ^(Tj 1 ^) governed by a 
potential force f{v,j) and a friction force proportional to 
dv/dx: 



d 2 v f da\dv 
K d^ + { CvS - Jd ^)di = f ^ Jd) 



da\ dv 



(20) 



(here x is the "time" and v is the "coordinate" of the 
particle). 

The velocity s of the domain is found from the condi- 
tion that the total change of energy E{x) (see Eq. (B3)) 
in the "period" L vanishes: 



\// c ( ! y)(^) 2 

\ dv V dx ) / \ dx 
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Hence, the velocity of a magneto-thermo-electric domain 
is 



ja/c v ~ j- 



1 fc s T c (1) fc B T c (1) 



c v T ( 



(i) 



10 cm/s 



for j - 10 B ^- 10 7 A/cm 2 , c v - 1 J/K cm 3 , T c (1) - 100 K. 

In the same manner as for Gunn domains in 
semiconductors and electric domains in superconductors 
and normal metals , the motion of the domain can be 
stopped by inhomogeneities inside the stack or at its ends. 
However, there is a possibility that the domain can move 
along the sample, periodically disappearing at one end of 
the sample and then reappearing at the other, so that the 
result is temporal non-linear electrical oscillations 29 with 
period uj ~ L/s. In the case under consideration these 
oscillations involve the magnetization direction 0, tem- 
perature T and electrical current j. Using the parameter 
values already introduced one gets uj = s/L ~ l-=-10 MHz 
for a magnetic stack of length L ~ 1 pm. 

As shown in Appendix B (see Eq. B4), the function 
v(x) satisfies the following equation with a high accuracy: 



1 



dx 



+ W(p,j d ) = E. 



(22) 



Here the constant E plays the role of "particle" energy, 
which determines the "period of motion", L(E,j), of a 
non-linear^oscillator. Its value is found from the condi- 
tion that L(E, j) is equal to the length L of the magnetic 
stack, that is the magneto-thermo-electric domain v{x) 
as a function of x has only one maximum and one min- 
imum along the length of the magnetic stack. In the 
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FIG. 6: Coordinate dependence of the temperature, T(x), 
and the magnetization- misorientation angle, Q(x), in a mag- 
netic stack containing a magneto-thermo-electric domain, cal- 
culated for R(Q) =R+-R- cos 6, R-/R+ = 0.2, D = 0.2; 

J/J c = 1.0265, J c = V c /R(tv), V c = y/RWQiTP)/^, 

L = [f c p{*)/^Y /2) ~L C . 



case L(E,j) = L/n (where n = 2,3,...), "multiple" do- 
mains are also possible in principle. However, they are 
unstable 14 and are therefore of no interest for us. 

The solution of Eq. (22), together with Eq. (17), de- 
fines the spacial distributions of temperature and mag- 
netization direction in a magnetic stack with a magneto- 
thermo-electric domain. Typical examples of such distri- 
butions are presented in Fig. 6 

Using Eq. (22) and the last equation in (19) one finds 
the set of equations 



L(E,j)=V2 
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Mj) y/E-W(T,j) 



= L, 
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The solutions E = Ed(V,L) and j = jd(V,L) of these 
equations are, respectively, the "energy" of the domain 
and the current-voltage characteristics of the magnetic 
stack containing a magneto-thermo-electric domain (the 
dynamic CVC). In Eqs. (23) and (24) the limits of inte- 
gration, i^ m i n and f max , i.e., the minimum and maximum 
temperature of the domain, are obtained as the real roots 
of the equation W(y, j) = E. 

Let us consider the most pronounced case, for which 
L S> L c . In this limit, as "multiple" domains are unsta- 
ble, one needs to find a solution of Eq. (22) that has a 
period much larger than L c . Since the period L(E,j) 
of the non-linear oscillator [see Eq. (23)] diverges loga- 
rithmically as E — >min{ Wi$ (j)} (where W\^{j) are the 
maxima and minima of the "potential" energy, see Fig. 4), 
this condition is satisfied if the "energy" E = E(L,j) of 



the domain differs from min{Wi, W3} by an exponen- 
tially small amount. From here it follows that variations 
of the current inside a narrow interval around jo , defined 
by the relation 

W 1 (j ) = W 3 (j ), 

drastically changes the form of the magneto-thermo- 
electric domain and hence the current voltage character- 
istics. Solving the set of equations Eq. (23) and Eq. (24) 
in the interval \j — jo\ "C jo (where W\{j) « Ws(j) ~ 
W(jo) ) and taking into account that the maximal, u max , 
and minimal, v m i n , temperatures of the domain are close 
to T\ and T3, respectively, one find the following implicit 
form of the dynamic CVC: 
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Here £ = V/L, L = [f T {Tr')/ n} (1/2) - L c and the 
constants Jjjj are of the order of j while /V = df /dT. 

Differentiating both sides of Eq. (25) with respect to £ 
one finds that the dynamic CVC have vertical tangents 
at the two points (Ji,£i) and (j2,£2), where 



and 
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j (l + A(7T)) 

p(7r)j {l-A(7r)[l-lnA(7r)} 



while 



A(6) 



L p (tt) - p(0) 



L 



p(Q) 



(27) 



(28) 



It follows from the above equations that the values of 
the currents ji and j'2 differ from jo by an amount ~ 
jo(L c /L). 

One can see from Eq. (25) and Eq. (28) that for all 
values of £ = V/L in the interval £ 1 < £ < £ 2 , ex- 
cept for a small region near the ends of the interval, the 
current j coincides with j to an accuracy that is expo- 
nential in the parameter L/ L c 3> 1. A change of the bias 
voltage in this interval does not change the current but 
it does change the length of the magneto-thermo-electric 
domain, that is the length of the higher resistive part of 
it providing the needed voltage drop across the stack at 
the fixed current value J — Jq. Therefore, a magnetic 
stack with a magneto-thermo-electric domain inside can 
work as a high-quality current stabilizer. 

Near the points £\ and £2 there is a sharp transition 
from the nearly horizontal segment of the dynamic CVC 
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to the rising segments of the CVC corresponding to a 
homogeneous state of the magnetic stack (see Fig. 5). As 
a result, there are hysteresis loops in the current-voltage 
characteristics of the stack of a large enough length L. 

The fact, mentioned above, that the current is nearly 
independent of the bias voltage in the interval [Vi , V2] 
(here Vi t 2 = L£\.2 comes about because in this region the 
domain structure can be written to exponential accuracy 
(i.e., with an error oc exp(— L/L c )) as 



Td(x) = 7?(a; + x\) + d(x 2 — x) — T n 



(29) 



where the function $(x) is a domain-wall type solution of 
Eq. (20) at j — jo and L — > 00, the asymptotic behaviour 
of which is 



lim = §{x) -> T n 
lim = <&(x) -> T m 



(30) 



Here x\,2 are the points of deflections of the curve 
dTd(x)/x so that Ljj = x 2 — x\ is approximately the 
length of the "hot" section of a trapezoidal MTED hav- 
ing the maximal temperature T max and Lj = L — L is 
the length of its "cold" section having the minimal tem- 
perature T m in (see Fig. 6) 



IV. CONCLUSION 

We have shown that Joule heating of the magnetic 
stack sketched in Fig. 1 by a current flowing in the plane 
of the layers may result in an instability of an initially 
homogeneous distribution of the magnetization of the 
stack if the length L of the stack in the direction of 
the current flow is longer than some critical length L c . 
This instability results in the spontaneous appearance of 
moving domains of magnetization direction in layer 2 of 
Fig. 1, 0(x — st), temperature, T(x — si) and electric 
field, £ (x — st). For the case L > L c the length of the 
domain is of the order of L c . 

If the length of the stack greatly exceeds the critical 
length, L 3> L c , the stack is spontaneously divided into 
two regions, where in one region the magnetization di- 
rections in layer 1 and layer 2 of Fig. 1 are parallel to 
each other, while in the other region they are antiparallel. 
The length of the region with antiparallel magnetization 
orientations (that is the length of the domain Lj) is con- 
trolled by the bias voltage in the interval L c < Lj < L 
(see Fig. 6). In this case the CVC of a stack containing 
such a domain has a plateau with hysteresis loops at the 
ends, as shown in Fig. 5. Therefore, the stack can work 
as a current stabilizer since the current flowing through 
it has a fixed value Jo to within an exponentially small 
error oc exp(— L/L c ); a change of bias voltage only results 
in a change of the domain length to provide the needed 
voltage drop over the stack. 

For a realistic experimental situation the value of the 
parameter r, which through Eq. (11) determines the de- 



pendence of the stack resistance on the magnetization- 
misorientation angle 0, can be estimated to be of order 
0.1 -T- 0.3, while the Curie temperature Tc of the spacer 
layer 1 can be ~ 100 K. Using these values and a re- 
sistivity of p(Q = 0) ~ 10 fxQcxn, a current density of 
j ~ 10 6 10 7 A/cm 2 and a specific heat of c v ~ l J/K 
cm 3 one finds that the critical length is L c ~ 10 /im. 
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Appendix A: Instability of a magnetization 
distribution that is spatially homogeneous in the 
plane of the magnetic layers 

In order to investigate the stability of the spatially ho- 
mogeneous solution (14) for temperature, T = To, cur- 
rent density, j = jo, and magnetization misalignment 
angle, = ©o, against spatial fluctuations we express 
these quantities as sums of two terms, 

T =T (V)+T 1 (x,t); Q = Q (V) + 0i(x, t) 
j =io+ji(*), (Al) 

where Ti, Qi and ji each is a small correction. Inserting 
Eq.(Al) into Eqs. (13) and (2) and using the Fourier 
expansion 

+00 



Ti(t,x) = Tl n >(t)exp{ik n x} 



n— — 00 
+00 



6i(t, x) — Qi(i)exp{ik n x}; h 



27m 



n— — 00 



one finds that 

T[ n \t) = DoTj^sineo/eo/ie^) 
while the equation for the Fourier harmonics of the angle 



de[ 0) 
Cv dt 

if n = 0, and 

c v dO^ 
K dt 



■{■ 



JcPM 



D T?> I (sine/0)' 
p(Q) (1 -7Jsin9/e) } 



(¥)-M- 



e=e, 



e^ 0) (A2) 



/t£> T c (1) |(sine/e)'| 
p(e)L- 1 (e)(l-7J sin 9/6)1 } e[ n) ; (A3) 

L J J W-Hn 
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if n + 0. 
From this result one sees that if 

[ P - 1 (e)(i-^sine/e)] e=eo <o 

[that is if dJ/dV < 0, see Eq. (9)] the uniform magnetiza- 
tion along the stack looses its stability if the stack length 
exceeds some critical value. Setting the right-hand side 
of Eq. (A3) equal to zero one gets the result (15) for the 
critical length L c . 



Appendix B: Domain structure 

By multiplying both sides of Eq. (20) by dv/dx one 
finds that the quantity 



E(x) 



(Bl) 



plays the role of the total "energy" of a fictitious par- 
ticle, the first term giving the "kinetic energy" and the 



"potential energy" being defined as 



W(T,j) 



f(T ,j)dT 



(B2) 



T2U) 



(the dependence of W(T, j) on T is shown in Fig. 4). The 
change of energy with "time" x is caused by the action of 
the "friction" force: 



dE 
dx 



(B3) 



As one sees from Eq. (20), the ratio of the "friction 
force" to the "incrtial" term is of the order 



OtL r 



rpK 



(B4) 



where sf is the Fermi energy. From here it follows that 
the function v(x) satisfies an "energy" conservation law 
E(x) = E = const (see Eq. (Bl)) to within an error 

~ fc B T c (1) /£F < 1, which results in Eq. (22). 
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